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Abstract 

We consider the region of small transverse momenta in the production of 
high-mass systems in hadronic collisions. By using the current knowledge on 
the infrared behaviour of tree-level and one-loop QCD amplitudes at (9(ag), 
we analytically compute the general form of the logarithmically-enhanced 
contributions up to next-to-next-to-leading logarithmic accuracy. By comparing 
the results with g^-resummation formulae we extract the coefficients that 
control the resummation of the large logarithmic contributions for both quark 
and gluon channels. Our results show that within the conventional resummation 
formalism the Sudakov form factor is actually process-dependent. 
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1 Introduction 



The transverse-momentum distribution of systems with high invariant mass produced in 
high-energy hadron colhsions is important for QCD studies and for physics studies beyond 
the Standard Model (see, e.g., Refs. ^-^). 

We consider the inclusive hard-scattering process 

hi{pi) + h2{p2)^F{Q\q^T-(P)+X , (1) 

where the final-state system F is produced by the collision of the two hadrons hi and 
/i2 with momenta pi and p2, respectively. The final state F is a generic system of non- 
strongly interacting particles, such as one or more vector bosons [j* ,W, Z, . . .), Higgs 
particles {H) and so forth. We denote by y/s the center-of-mass energy of the colliding 
hadrons {s = {pi +^2)^ — '^PiP2)i and by and the invariant mass and total transverse 
momentum of the system F, respectively. The additional variable in (0) denotes the 
possible dependence on the kinematics of the final state particles in F (such as rapidities, 
individual transverse momenta and so forth). 

We assume that at the parton level the system F is produced with vanishing qx (i.e. 
with no accompanying final-state radiation) in the leading-order (LO) approximation. Since 
F is colourless, the LO partonic subprocess is either qjqji annihilation, as in the case of 
7*, W and Z production, or gg fusion, as in the case of the production of a Higgs boson H . 

When the transverse momentum of the produced system g|, is of the order of its invariant 
mass the fixed order calculation is reliable^. In the region q^ <^ large logarithmic 
corrections of the form a^/q^log^""'^ /q^ appear, which spoil the convergence of fixed- 
order perturbative calculations. The logarithmically-enhanced terms have to be evaluated 
at higher perturbative orders, and possibly resummed to all orders in the QCD coupling 
constant as- The all-order resummation formalism was developed in the eighties 0- 



The structure of the resummed distribution is given in terms of a transverse-momentum 
form factor and of process-dependent contributions. 

The coefficients that control the resummation of the large logarithmic contributions 
for a given process in (0) can be computed at a given order if an analytic calculation at 
large at the same order exists. At first order in as the structure of the large logarithmic 
contributions is known to be universal and depends only on the channel in which the system 
is produced in the LO approximation. At second relative order in as, only a few analytical 
calculations are available, like the pioneering one for lepton-pair Drell-Yan production, 
performed by Ellis, Martinelli and Petronzio in Ref. |]15[. Using the results of Ref. |T5| 



Davies and Stirling (see also [|l6l) were able to obtain the complete structure of the 
(9(ag) logarithmic corrections for that process. 

The analysis performed by Davies and Stirling is by far non trivial because it requires the 
integration of the analytic qT distribution in the small limit. Moreover the calculation 
cannot tell anything about the dependence of these coefficients on the particular process 
in (|1|) and should in principle be repeated for each process. 

^It is assumed that all other dimensionful invariants are of the same order Q^. 
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In this paper we address this problem with a completely independent and general 
method. Our basic observation is that the large logarithmic corrections are of infrared 
(soft and coUinear) nature, and thus their form can be predicted once and for all in a 
general (process independent) manner. 

The structure of the logarithmically-enhanced contributions at 0(ag) is controlled by 
the infrared limit of the relevant QCD amplitudes at the same order. The infrared behaviour 



of QCD amplitudes at 0{as) is known since long time ^Tj. Recently, soft and collinear 
singularities arising in tree-level |T8|, |19[ and one-loop |21], QCD amplitudes at 

(9(ag) have been extensively studied and the corresponding kernels have been computed 
|T^-|2^. By using this knowledge, and exploiting the relatively simple kinematics of the 
process we will construct general approximations of the relevant QCD matrix elements 
that are able to control all singular regions corresponding to g-r ~^ avoiding double 
counting. By using these approximations we will compute the general structure of the 
(9(ag)-logarithmically-enhanced contributions both for qq- and for 5f(yf- initiated processes. 

The results provide an important check of the validity of the resummation formalism 
and allow to extract the general form of the resummation coefficients. In particular in the 



quark channel we can confirm the results of Ref. [jT4| in the case of Drell-Yan and in the 
gluon channel we can give the coefficients in the important case of Higgs boson production 
through gluon-gluon fusion. 

The universality of our method relies on the fact that the infrared factorization formulae 
we use depend only on the channel [qq or gg) in which the system F is produced at LO 
and not on the details of F. 



Our main results were anticipated in a short letter ||2g]. This paper is organized as 
follows. In section 2 we review the framework of the resummation formalism and present 
the strategy for the calculation. In section 3 we perform the calculation explicitly for the 
0{as) corrections and extract the first order coefficients. Section 4 and 5 are devoted to 
the calculation at 0{ag) for the quark and the gluon channel and constitute the main part 
of this work. Finally in section 6 we present our final results and discussion. 



2 Resummation formula 



The transverse momentum distribution for the process in Eq. (|I|) can be written as: 

dap 
dQ"^ dq^ d(j) 

Both terms on the right-hand side are obtained as convolutions of partonic cross sections 
and the parton distributions fa/h^x^Q"^) (a = qf,qf,g is the parton label) of the colliding 
hadronsQ. 

The partonic cross section that enters in the resummed part (the first term on the 
right-hand side) contains all the logarithmically-enhanced contributions a'^/q^log"^ Q"^ /q^. 

^Throughout the paper we always use parton densities as defined in the MS factorization scheme and 
as(g^) is the QCD running couphng in the MS renormahzation scheme. 



dar 



dQ"^ dq^ 



+ 



dar 



dQ"^ dqi. 



(2) 



. fin. 
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Thus, this part has to be evaluated by resumming the logarithmic terms to all orders in 
perturbation theory. On the contrary, the partonic cross section in the second term on the 
right-hand side is finite (or at least integrable) order-by-order in perturbation theory when 
— ^ 0. It can thus be computed by truncating the perturbative expansion at a given fixed 
order in as- 

Since in the following we are interested in the small-gr limit we will be concerned only 
with the first term in Eq. (|^). The resummed component is0 



dap 
dQ"^ dq^d 



/ dxi j dx2 db - Jo{bqT) fa/hi{xi,bl/b'^) fb/h2{x2,bl/b^ 
„u Jo Jo Jo ^ 



■ sPF„^(xiX2s;(5,6, 



(3) 



The Bessel function Joi^bgr) and the coefficient 60 = 2e~'''^ (7^ = 0.5772. . . is the Euler 
number) have a kinematical origin. To correctly take into account the kinematics constraint 
of transverse-momentum conservation, the resummation procedure has to be carried out in 
the impact-parameter 6-space. The resummed coefficient is 



= E 

c 



' dz, f dz, C^Mbl/b'), z,) C^,{as{bl/b'), z,) 5{Q' 
Jo 



Z1Z2S) 



(4) 



where da^^*^^ /d(j) corresponds to the leading order cross section for the production of the 
large invariant mass system F in the cc channel, with c representing either a quark g or a 
gluon g. The resummation of the large logarithmic corrections is achieved by exponentia- 
tion, that is by showing that the Sudakov form factor can be expressed as 



SciQ, 6) = exp I - [ 
[ J^l 



dq 



2 r 



A,(as(g')) log^ + i?,(«s(g')) 



(5) 



The functions Ac{as)^ Bc{as), as well as the coefficient functions Cafe(as, z) in Eqs. (§J^) 
are free of large logarithmic corrections and have perturbative expansions in as as 

00 



n=l 

00 



n=l 



00 

CUas,z) = S^bSil -z) + J2{^ycj-\z) . (8) 



n=l 



The coefficients of the perturbative expansions yir, grander 

(z) are the key of the 

resummation procedure since their knowledge allows to perform the resummation to a given 

■^This expression can be generalized to include the dependence on the renormalization and factorization 
scales fiR and fip, respectively (see e.g. Ref. p7|). 
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Logarithmic order: A^^^ leads to the resummation of leading logarithmic (LL) contributions, 
{y4(2),i?(^),C(^)} give the next-to-leading logarithmic (NIL) terms, {A^^\ B^^\C'^'^^} give 
the next-to-next-to-leading logarithmic (NNLL) terms, and so forthQ. The coefficient func- 
tions Cj^^ (z) depend on the process, as it has been con&med by calculations of C^^^ (z) for 
several processes. The Sudakov form factor Sc{Q, b) that enters Eq. (^) is often supposed to 
be universal. However, as we will show, this is not the case, and anticipating our results we 
label all process-dependent coefficients by the upper index F. The coefficients A^^\ B^^\ 
are universal and are known both for the quark [jlO[ and for the gluon |jl3|] form factors 



^(1) 


= 2Cf , 


41) 


= 1Ca 




= — 3Cf , 


bw = 


-2/?o, 




= 2CfK , 




2CaK 



(9) 



where 



and 



/?o = yCA-^%r^ (10) 



The NNLL coefficient B^^) was computed by Davies and Stirling W% for the case of Drell- 
Yan (DY): 

5f = Cl{.'-\- 12C(3)) + Cf Ca - ^ + 6C(3)) + Cf Tn (j - ^tt^) 

(12) 

where C(n) is the Riemann (^-function (C(3) = 1.202 . . . ). It is also worth noticing that, even 

(3) 

though there is no analytical result available for it, the coefficient Ag^g has been extracted 



numerically with a very good precision in Ref. p9|. 



As anticipated in the introduction, a direct way to obtain the coefficients in Eqs. 
at a given order involves the computation of the differential cross section da / dq^dQ"^ d(j) 
at small qt at the same order. A comparison with the power expansion in as of the 
resummed result in Eq. (^ allows to extract the coefficients that control the resummation 
of the large logarithmic terms. However, it has been shown by Davies and Stirling that is 
it more convenient to take z = Q'^/s moments^] of the differential cross section defining the 
dimensionless quantity 

^^^^^Jo '^'''''d^iJd^d^dQ^- ^^^^ 

Notice that in the definition of S the cross section has been normalized with respect to the 
lowest order partonic contribution daQ/d(f) and multiplied by to cancel its l/q^ singular 

''in a different classification tlie coefficient C'^^ enters only at NNLL [2^ . 

^Here we follow Ref. in tlic miconvcntional definition of tlic moments: f{N) — dzz^ f{z). 
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behaviour in the hmit Qt ^ 0. The upper hmit of integration z = 1 — 2^^/(5(^/1 + C[t/Q'^~ 
Qt/Q) ~ 1 ~ '^Qt/Q has been approximated to a first order expansion in qx/Q and cor- 
responds to the kinematics for the emission of soft particles (i.e, when the center of mass 
energy s is just enough to produce the system with invariant mass Q and transverse momen- 
tum qT). Working with moments allows to avoid complicated convolution integrals implicit 
in d^) and makes possible to factorize the parton densities from the partonic contribution 
to the cross section. In this way, the corresponding expression from the resummed formula 
(I) reads 

E(iV) = Y,u,nAN,^^l)f,,hAN,^^l)T.,m (14) 

where 

poo 2 



a,b 







■ exp 



2 



dq 

62/fe2 q^ 



2 

2 



A,(as(g^)) log^ + i?f(as(g^)) 



i2 



(7ai + 7bj)(A^,as(r)) 



(15) 



and an ordered exponential is understood. Notice that the appearance of an extra term 
involving the anomalous dimensions 7af) in the exponential in (|15]) is due to the evolution of 
the parton densities from the original scale h'^/h'^ in (|^) to the arbitrary factorization scale 
/ii? at which they are now evaluated. 

In order to extract the resummation coefficients, we can directly study the partonic 
contribution Sjj. Furthermore, since we want to perform a calculation of Ej^ to 0{a^ 
and our main interest is the second order coefficient B'^'^\ it is clear that only the diago- 
nal contribution Ecg can give the desired information. Each possible "flavour changing" 
contribution in Eq. (|15D would add at least one extra power of as in the perturbative 
expansion. 'Non- diagonal' contributions to Sjj, which can be evaluated in a simpler way, 
might be used to check the structure and consistency of the resummation framework at a 
given perturbative order but do not provide any additional information on the coefficients. 

In order to have transverse momentum 7^ at least one gluon has to be emitted and, 
therefore, the perturbative expansion of Sec begins at 0{a^) 

Sc,(iV) = ^ (N) + (^) ' Eg) (iV) + . . . . (16) 



From the expansion of the resummed formula (|I5|) it is possible to obtain the expression 
for the first two coefficients in (|16]) a^ 

E«(iV) = A« log ^ + + 27L^)(iV) (17) 



®For the sake of simplicity in the presentation, and unless otherwise stated, we fix the factorization and 
renormalization scales to = /i^ = Q^. 
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and 



+ log" \ 

+ log ? + {B^^ + 27£)(iV)) - {Bi'^ + 2l^HN)y 

Qrp L 



- 2/5oC£)^(iV) + 2 5^ fc«^(iV)7];)(iV)" . 



The computation of (A^) can provide information on the first order coefficients A^P 
(the logarithmic term in ([T7|) ) and B^^ (the constant term in (|l3)) as well as on the one- 
loop anomalous dimensions 'ycl\N) ( the A^— dependent term in (p^TD H- ) ■ In the same way, 
the coefficients A^c^ and B^^^^ can be extracted from the second order result (ITSf). At this 



order, also the coefficient functions Clj^^{N) contribute to the logarithmic and constant 
terms and therefore should be known in order to be able to proceed with the extraction of 
Ac and Be . Fortunately, there is another related quantity which allows to obtain the 
coefficient functions C^y (A^) from a first order calculation. This is the gr-integrated cross 
section 



2 ^ic ■ 



(19) 

/o It 

When p\. <^ the perturbative expansion to 0(q;s) reads (neglecting again terms that 
vanish when pt — ^ 0) 



-y^^ log^ % - {B^^ + 2lil\N)) log ^ + 2C(r (AT)' 



2--C -o ^2 V c /cc V JJ -6^2^ 



as 

27T 



Q 



Pt 



i ^ c. 



(20) 



The integration over q^ adds one power in the logarithm, with the coefficient functions 
C^j^^{N) appearing now in the constant term. It is worth noticing that at variance with 
the calculation of S the configuration with = now contributes to Eq. ([TP|). 

In the quark channel (c = g), for the sake of simplicity and in order to compare directly 
with the calculation performed in |T^, we will concentrate on the non-singlet contribution 
to the cross section defined by 



NS _ \^ / _ 

ff 



(21) 



^Notice that all moments but one can actually be extracted. The remaining one can be obtained by 
imposing quark number and momentum conservation rules. 
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The second order expansion for S^^(A^) in terms of the resummation coefficients reads 
hke the one in Eq. (0) but without the 'singlet' contributions involving Ylj^c ^^"^ with 
the corresponding non-singlet anomalous dimension. In the following the label NS will be 
always understood in T,qq. 



3 The calculation at 0{as) 

The calculation at 0{as) is not difficult and the results are rather well known. Nevertheless, 
we will give in this section the details on the computation as a way to present the main 
ideas of the method developed to obtain the resummation coefficients. 

At this order only one extra gluon of momentum k can be radiated and the kinematics 
for the process cc g + F is (see Fig. |l|) 

Pi+P2^k + q. (22) 




Figure 1: 0{as) contribution to the process (1) 

We denote the corresponding matrix element by A^c^L^f (Pi)P2, 0) and the usual 
invariants are defined as 

s = {pi+P2f u = {p2-kf t={pi~kf z = Q^/s. (23) 

The differential cross section can be written as 

d^^^9F _ f \M^S-.,FiPuP2,k,^)\'iAnyq-^^ du (1 \ 
dqldQHcl>~ J 8.(27r)2 1(1-6) 7' 

where the two roots of the equation {pi + P2 — q)"^ = ^ are given by 



Iz Iz 

(25) 

In order to regularize both ultraviolet and infrared divergences we work in the conventional 
dimensional regularization scheme (CDR) with 4 — 2e space-time dimensions, considering 



7 



two helicity states for massless quarks and 2 — 2e helicity states for gluons. The lowest-order 
cross section (at = 0) needed to construct S in Eq. (|T3|) is given by 



rf0= 2s ' ^26) 

in terms of the Born matrix element |A^^5^f (PijPa, 4>) 

As has been stated, we want to obtain S^g'' by using our knowledge on the behaviour 
of QCD matrix elements in the soft and collinear regions at (9(ag). The starting point 
is the observation that, when is small, the additional gluon is constrained to be either 
collinear to one of the incoming partons or soft. Thus there are three singular regions of 
M^^~_^p (pi,P2, k, 0) in the qr ^ limit: 



• first collinear region: pik 

• second collinear region: — >■ 

• soft region: /c — > 0. 



It is clear that, since q'^ is small but does not vanish, these regions do not produce any real 
singularity, i.e. poles in e, but are responsible for the appearance of the logarithmically- 
enhanced contributions. When pi/c — the matrix element squared factorizes as follows: 



\M 



(0) 

cc^gF 



(,Pl,P2,k, 



zipik 



2e 



where 



Pqq{z^ e) — Cp 
ICa 



Pa9{z,e 



1 + 
1 - z 

z 

1 - z 



+ 



e(l 
1 - z 



z 



+ z{l-z) 



(27) 

(28) 
(29) 



are the e-dependent real Altarelli-Parisi (AP) kernels in the CDR scheme. In the left hand 
side of Eq. ( ^7]) the matrix element squared is obtained replacing the two collinear partons 
c and (7 by a parton c with momentum Zipi. 

Notice that in the gluon channel there are additional spin-correlated contributions and 
Eq. (p7|) is strictly valid only after azimuthal integration. Since here and in the following 
we will always be interested in azimuthal integrated quantities, Eq. (|2^) can be safely used 
also in the gluon channel. 

In the limit p2k ^ the singular behaviour is instead 



\M 



(0) 
cc^gF 



iPuP2,k,(f)) 



Z2P2k 



-Pcc{z2,e)\M 



^FiPU^2P2, 



(30) 



Let us now consider the limit in which the gluon becomes soft. As it is well known soft- 
factorization formulae usually involve colour correlations, that make colour and kinematics 
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entangled. In general colour correlations relate each pair of hard momentum partons in the 
Born matrix element. In this case the hard momentum partons are only two and colour 
conservation can be exploited to obtain: 



where 



is the usual eikonal factor and we have defined 

Cq = CF Cg = CA- (33) 

In Eq. (^) colour correlations are absent and factorization is exact. This feature will 
persist also at C(ag). 

In each of the singular regions discussed above, Eqs. (^7|), (|30|) and (|3l|) provide an 
approximation of the exact matrix element that can be used to compute the cross section 
in the small qt limit. In principle it might be possible to split the phase space integration 
in regions where only soft or collinear configurations can arise, and use in each region the 
corresponding approximation. Unfortunately, such method probes to be very difficult to be 
extended to (9(a;g), where the pattern of singular configurations is much more complicated. 
Thus our strategy is to unify the factorization formulae in order to obtain an approximation 
that it is valid in the full phase space. 

As can be easily checked, if we identify the momentum fractions zi and Z2 with z, the 
collinear factorization formulae in Eqs. (|27|j30| ) contain the correct soft limit in Eq. (PH). 
Therefore, the unification of soft and collinear limits is rather simple: the usual collinear 
factorization formula already contains both. Strictly speaking, one can use the symmetry 
in the initial states in order to perform the integration in Eq. (p^ only over half of the 
phase space (i.e. by taking for instance only u = Umax) and multiplying the result by two. 
In this way only one possible collinear configuration can occur and Eq. (p7|) provides the 
needed approximation for the matrix element. 

At this order it is even possible to write down a general factorization formula for the 
three configurations that shows explicitly the I/qt singularity of the matrix element squared 

as 

\M^S-.,APuP2,k,<P) P ^^^^Ae(.,e)|A<2U(0) \\ (34) 

where we have used Lorentz invariance in order to write \M^^^_^p{(f))\'^ only as a function of 
the final state kinematics. We can now use this formula to compute the small qt behaviour 
of Scc(iV) in a completely process independent manner. In fact the process dependence, 
given by the Born matrix element, is completely factored out and cancels in E. By replacing 
Eq. (|3^ ) in Eq. ( ^4]) and using the definition of S we obtain, keeping for future use its e 
dependence: 



S«(iV,e) = -— -(^) / dzz 



1 fAnfi^Yf'-'^-/'^^ ^ 2{1- z)PUz,e) 



r(l -e)\ ) J, v/(l-^)2-4zg2/Q 



2 



1 Avr/i^ 



C,T,^N,e). (35) 
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Explicitly, setting e to 0, we have 

s(^,)(iV) = 2Cf log ^ - + 27;;) (iV) 

S« (iV) = log ^ - 2/3o + 27^) (iV) , (36) 

for the quark and gluon channels. Comparing to Eq. (|1^ we see that A^c"* = 2Cc is 
the coefficient of the leading 1/(1 — 2) singularity in the AP splitting functions whereas 
is given by the coefficient of the delta function in the regularized AP kernels 

1? = li'^ = Po . (37) 



Finally, in order to obtain the coefficient C^|,\ we have to evaluate the integrals in 
Eq. (|19D and compare to the results from Eqs. (|20D. As far as the diagonal contribution is 
concerned, one has to take into account also the one-loop correction to the lowest order cross 
section, a contribution formally proportional to S{q^). The interference between the one- 
loop renormalized amplitude with the lowest order one depends of course on the process. 
Nevertheless, its singular structure is universal and allows to write in general [BD 



27r V J r(l - 2e) V e 

The finite part depends (in general) on the kinematics of the final state non-coloured 
particles and on the particular process in the class (|l|) we want to consider. In the case of 
Drell-Yan we have [^]: 



= CF{-8 + ln'] , (39) 



whereas for Higgs production in the rritop 00 limit the finite contribution is |^ 

2 
3 



Af = 5Ca + -Catc^ -3Cf = 11 + 2n\ (40) 



The diagonal term in Eq. ([T9|) can be evaluated integrating Eq. (p5|), from to p^, keep- 
ing into account the contribution in Eq. (|38| ) and subtracting the following factorization 
counterterm in the MS scheme: 

As for the non-diagonal contribution, one needs T,ic{N), that can be computed, analogously 
to Eq. (ID as 



r(l -e)\ ql ) v/(l-^)2-4zg^/Q2 



1 / 47r/i 



r(l -e)\ql 



dzz^'P^^iz^e) (42) 
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where the functions Pci{z,e) are the non-diagonal AP sphtting kernels 



1 + (1 



T, 



R 



Z 

2z{l 



ez 



and the absence of singularities as z 



1-e 

1 has been exploited to set qr 



The factorization counterterm to be subtracted in this case is 



R 



(FCT) 



(iV) = - 



1 rfi 



rfl -2e) 



(N) 



(43) 
(44) 

in the integral. 

(45) 



Comparing the total results to Eqs. 



we obtain for C, 



(1). 



-KiU) + i^Hi - ( a -g- + 2-4: w 



TT 



ab 
2 



(46) 



where Pab{z) represent the 0{e) term in the AP Pab{z,e) splitting kernels in Eqs. 
|3|, and are given by: 



, 29 



Pg,{z) 
P'Jz) 



(l-z) 



-Cf(1 
-Cfz 
-2Trz{1 
0. 



(47) 



As can be observed, the coefficient function contains both a hard process dependent contri- 
bution (proportional to (0)) originated in the one-loop correction as well as a 'residual' 
collinear contribution proportional the e part of the splitting functions which has origin 
in the particularities of the MS scheme (see Eq. (^iD), where only the e = (and not the 
full) component of the splitting functions is factorized. The general expression in Eq. (^) 
reproduces correctly the coefficient C^J^^ computed for Drell-Yan |T^, Higgs production in 
the rritop oo limit |^ 0, 77 and ZZ |^ production. 



Summarizing the 0{as) results, the coefficients A^c'^ and Bc'^' are fully determined by 
the universal properties of soft and collinear emission. The function C^^'* depends instead 
on the process through the one-loop corrections to the LO matrix element. 



?(i) 



4 The calculation at 0{al): the quark channel 

At 0{al) T,qq{N) receives two contributions: 

• Real emission of two partons recoiling against the final state system F^Q"^, qT,4>)l 

• Virtual corrections to single-gluon emission. 

In the following we compute these contributions in turn. 
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4.1 Real corrections 



The computation of the double real corrections to T,qg{N) represents the most involved part 
of the complete calculation. The difficulties arise both from the fact that the additional 
parton in the final state implies three more phase space integrals, and from the appearence 
of many more singular configurations that contribute to the limit qx 0. 



The kinematics for the double real emission process cc ^ i + j + F is (see Fig. 

Pi+P2^ h + h + q, 



(4J 




^ F 



Figure 2: 0{ag) contribution from double real emission. 

and the corresponding matrix element is denoted by f{pi,P2, ki, k2, 

usual invariants are defined as 

s = {pi+ P2f t={pi- qf u = {p2- qf S2 = {ki + ^2)^ , 



The 



(49) 



and fulfill the following relations 

S + t + U = Q'^ + S2 



2 Ut - S2Q'^ 

Qt = 



(50) 



In terms of these invariants, the real contribution to the cross section at fixed q^ is 
given by 



dq'^dQ^dcf) 
where dQ is 



|-M2U-^(pi,P2,fci,A;2,0)p 



{s2q'T)-' 



du dVt 



2s 



(47r)4-2T(l - 2e) Q'^-u2'k 



(51) 



dVL = sin^^^'d^2 smel^'^'ddi 



(52) 



with the angles defined in the frame where the partons corresponding to momentum ki and 
k2 are back-to-back ITH . 
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We see from Eq. (51) that the first step of the calculation involves the integration over 
the two angles. The integrals needed here are typical of heavy quark production at NLO 
and most of them can be found in Refs. |38| . The results of the angular integrals contain 



poles up to 1 /e while terms that develop an extra additional singularity as S2 — have to 
be computed up to C(e). 

The second step is the integration over u (or 52)- The integration limits are given by 
the two roots Umax and Umin in Equation (125|) . At this point, it is convenient to define the 
'symmetric' value for which u = t 



Uo = Q'- sjs{ql + g2) . (53) 
This value of u corresponds also to the maximum of S2 



and, in the CM frame of and p2, to the configuration where = 0. The singularity in 
S2 is made manifest by use of the identity 



s^'-' = —6{S2) 1 - e log A + -e^ log^ A + —— - e ^ + 0{e') , (55) 
e V 2 y (s2)a+ V S2 J A+ 

with the distributions defined as: 



ds2j^= / — (/(32)-/(0)) (56) 

[S2)a+ Jo S2 



r<i../M(i^) = r rf.. - /(O)) . (57) 

^0 V ■52 Jo S2 

In order to obtain S(A^) one finally has to integrate over z, keeping only the terms 
that do not vanish in the small-gr limit. At the beginning we consider only the N = 
moment^. We will later show how to perform the calculation for general A^, once one 
moment is known, in a simpler way. Notice that after implementing the regularization of 
the S2 = singularities using Eq. (^5l), the last two integrals can be performed directly 
in four-dimensions, since the small transverse momentum qt acts as a regulator of other 
possible singularities. 

(2) 

The double real contributions to (the non-singlet part of) S^^^(A^) fall into three classes, 
according to the possible different final states: 

• q + q^q + q + F 

• q + q-^q + q + F 

^Notice that the calculation of a single moment is enough to obtain the resummation coefficients A^^^^ 
and ^(2). 
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• q + q^g + g + F . 

Notice that the q + q^q + q + F is needed to form the non-singlet combination. 

As we did at 0{as), to study the small gr behaviour of S we will rely on the structure 
of soft and coUinear singularities of the corresponding QCD matrix element. In principle 
there are, of course, configurations where the two final state partons are hard and emitted 
back-to-back with small total transverse momentum. Nevertheless, these configurations do 
not produce any singularities when g^- — > and thus may be neglected. Finally, notice that 
we consider only double singularities, i.e. configurations where two extra partons are either 
collinear or soft, without caring about single singularities. Configurations with only one 
collinear or soft parton (and the other hard) do not contribute to S*^^) since the system F 
is not emitted with small qx in such case. 



4.1.1 Contribution from qq and qq emission 



For the q + q-^q + q + F contribution we have three singular regions at 0{ai^) [T^ 



• first triple-coUinear region: kipi ~ k2Pi ~ A;iA;2 ^ 0; 

• second triple-coUinear region: kip2 ~ k2P2 ~ ^1^2 0; 

• double-soft region: ki,k2^0 . 

In the first region the singularity is controlled by the following collinear factorization for- 



mula 18, 24, 19 



\M'-^l,,FiPl,P2,h,k2A)\' ^ ^^''^'^"^^' p,^,,,.(,3)I-^S^f(^3Pi,P2,0)P, (58) 

where Pq^q^q^t^q^) is the splitting function that controls the collinear decay of an initial state 
quark of momentum pi into a final state quark- antiquark pair qiq2 of momenta ki and 
^2 and the 'off-shell' quark q^ that participates in the hard cross section.. The explicit 
expression of Pq^q^q^^q.^^) is obtained from the one of Pq^q2qz, the splitting function for the 
decay of a ('off-shell') quark into a final state quark- antiquark pair plus a quark, given in 
Eq. ( |A.1|) , with the following definitions 



S12 = S2 Si3 = -2piki S23 = -2pik2 

Xi = -Z1/Z3 X2 = -Z2/Z3 X3 = 1/Z3, (59) 

where zi and Z2 are the momentum fractions of qi and q2 {z^ = 1 — Zi — Z2). Notice that 
Eq. ( p9D corresponds to the following transformation: 



n ki, r2 k2, rg -pi , (60) 
applied to the expression in Eq. ( |A.1| ) to cross the 'off-shell' parton to the final state. 
A formula similar to Eq. (p8|) can be written in the second collinear region. 
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In the double-soft region the factorization formula is instead 



19 



where 

[hhrPiih + k2) Pjih + k2) 
The reader can easily check that by defining the momentum fractions in Eq. (^9D as0 

hP2 k2P2 

zi = Z2 = (63) 

P1P2 P1P2 



Eq. (pSf ) correctly keeps into account also the double-soft limit in Eq. (pll). Thus, at least 
outside the second collinear region, the factorization formula (|58| ) with the definitions (|6^) 
correctly gives the full singular behaviour in this channel. 

The strategy to perform the calculation is the following. We use Eq. (|8]) to approximate 
the matrix element in its region of validity and compute its contribution to Eq|'*^(0) by 
integrating only in half of the phase space, that is from uq to Umax- The remaining region, 
which is obtained by exchanging u ^ t, will give , due to the symmetry of the initial state, 
exactly the same contribution and it is taken into account by multiplying the computed 
result by 2. As it happens at leading order, the information on the process, embodied in 
the Born matrix element is completely factored out in the calculation and disappears in E. 
In fact the Born matrix element can be fully written in terms of the (fixed) kinematics of 
the final state particles \Mf^^pizsPuP2, (j))]^ = |-MJ°L^(0)|'. 

For the q + q—*q + q + F contribution, needed to form the non-singlet contribution in 
Eq. (pID, there are only two singular configurations: 



• first triple-coUinear region: kipi ~ k2Pi ~ kik2 0; 

• second triple-coUinear region: kiP2 ~ k2P2 ^ kik2 . 

For the first collinear region we can write: 

\Mf^^^^APl,P2,h,k2,<l>)\' ^ ^3P2, 0)P • (64) 

Here -Pg^giqaCgs) is now the splitting function which controls the collinear decay of an initial 
state quark into a final state qq pair. The explicit expression for Pq~*qiq2(q3) can be obtained 
from the expression of Pqiqzqa in Eq- ( A.l ) with the following definitions 



S12 = -2pik2 Si3 = -2piki S23 = S2 

Xi = l/Z3 X2 = ~Z2/Z3 X3 = -Zi/Z3, (65) 



^To parameterize the triple-collinear limit it is necessary to introduce an additional light-cone vector n. 
This definition corresponds to the choice n = p2. Notice that a similar definition can be adopted also at 
0{as) to reobtain Eq. in the small qt limit. 
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i.e., corresponding to the crossing transformation: 



ri -pi, r2 k2, ki, (66) 

and similarly for the second coUinear configuration (with pi <-> P2). There is a partial 
cancellation between the CfTr contribution to E^^ from Eqs. ( pSj ) and (|6lD, due to the 
non-singlet combination. Once this cancellation is carried out, the part corresponding 
to the production of 'non-identical' partons in the qq channel gives the following CpT^ 



(2) 



contribution to 



5-^«-(0,<)-5iog^|-5iog| + i + -.= 



(67) 



where 



= rir^z^) (irJ W ■ 

and the explicit expression of function jFqg(0,e), defined in Eq. ( pSf ) is 



^,,-(0,e) =21og^-3-e. (69) 

At the beginning of Eq. ( |67D we have isolated a divergent term which will be cancelled by 
a similar one appearing in the virtual contribution. 

The part corresponding to the production of 'identical' partons in the qq channel gives 
also a Cf{Cf — Ca/'^) contribution, which does not contain any logQ^/g|. term. Therefore, 
there is a great simplification in the calculation since can be set to zero just after 
performing the angular integrations. We find: 

sS\ti)(0) = Cf{Cf - ic^) ( - 6 + 2vr2 - 16C(3)) . (70) 

The calculation of the qq contribution can be performed with exactly the same strategy as 
for the qq channel]^. After the CfTr contribution has been cancelled with a similar one in 
the qq channel only a contribution proportional to Cf{Cf — Ca/'^) remains[^ 

SS(^,)(0) = -Cf{Cf - ]-Ca) f ^ - + 4C(3)) . (71) 



4.1.2 Contribution from gg emission 

(2) 

The calculation of the double-gluon emission correction to S]^ is more difficult, because it 
is not possible to keep into account all possible singular configurations by using only the 
triple-collinear splitting functions. We will divide the calculation in two parts, according 

factor 1/2 has been included to account for the two identical particles in the final state. 
^^The overall minus sign here is due to the fact that this quantity must be subtracted in order to construct 
the non-singlet combination. 
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to the corresponding colour factors. First we will consider the non-abelian, CfCa term, 
which turns out to be simpler, and finally the abelian, part. 



CfCa contribution 



For this colour structure there are three singular regions to be considered [|T9 



first triple-collinear region: kipi ~ ^2^1 ~ A;i/c2 ^ 0; 
second triple-collinear region: kip2 ~ ^2^2 ~ ^1^2 — ^ 0; 
double-soft region: fci, A;2 — ^ . 



We point out that, as discussed in Ref. |T9[, thanks to the coherence properties of soft- 
gluon radiation, the soft-collinear region does not give any contribution proportional to 
CfCa (see later). 

In the first coUinear region the singularity is controlled by the following factorization 
formula: 



(0) 



{pi,P2,ki,k2,(l))\ 



nab 



CFCAP^':Z.ijKlFi^^P^^P2, 0)P , (72) 



where P*-"^^^-* is the non-abelian part of the splitting function that controls the collinear 
decay of an initial state quark into a final state gluon pair. This function can be obtained 
from Eq. ( |A.?| ) with the replacement in Eq. (|59|) . 



A similar formula to Eq. (^) can be written in the second collinear region (by pi ^ p2 
exchange) . 

In the double-soft region the factorization formula is instead (see Eq. (A. 3) of Ref. ||19|| ): 

\Mf^^ggFiPl,P2,ki,k2,(l))\l,y, ~ {ATTH^'as)^CFCA 

■ (25i2(fci, k2) - Su{ki, k2) - S22{ki, k2)) \Mqq-.F m\ 

(73) 

where the non-abelian double-soft function reads 

(1 - e) piki pjk2 + Pik2 Pjki 



Sijiki, k2) 



{kik2Y Pi{ki + k2) Pj{ki + k2) 



{PiPjf 



2piki pjk2 Pik2 Pjki 



2 - 



Piki Pjk2 +Pik2 Pjki 
Pi{ki + k2) Pj{ki + k2] 



(74) 



+ 



PiPj 
2kik2 



+ 



Piki Pjk2 Pjki pik2 



4 + 



Pi{ki + ^2) Pj{ki + k2 



{piki Pjk2+Pik2 Pjki 
Piki Pjk2 Pik2 Pjki 



As it happens in the qq and qq channels, it turns out that by defining the momentum 
fractions of the gluons as in Eq. (63), the factorization formula in Eq. (^) correctly accounts 
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also for the double-soft configuration. Furthermore, we have verified that Eq. (^) does not 
introduce any additional spurious singularities in the other infrared configurations. Thus 
for this colour structure the situation is similar to the one in the qq and qq channels and we 
can follow the same strategy. We approximate the non-abelian part of the matrix element 
in the region from uq to Umax using Eq. ([721). We first perform the angular integrals and 
then, exploiting the pi ^ p2 symmetry, do the remaining u and z integrations only over 
half of the phase space, i.e. with u from uq to Umax- 

In order to perform the last two steps, that are considerably more complicated than in 
the case of qq emission, we developed Mathematica programs that are able to handle 
the cumbersome intermediate expressions in the small qx limit. 



The result is0 



-'qqiggnah) 



(0) = CfCa JC 



11 , Q' 

— + log — 



6 



+ log^ ^ 



13 , 2 



35 



TT 

3 



-^,,-(0,6) 



log^ + 4C(3) 



2 + — tt' 
18 



(75) 



in agreement with Ref. ||0|. The first line of Eq. ( ffSD comes from the singular 6{s2) terms 
and will be exactly cancelled by a contribution appearing in the virtual correction. 



Cp contribution 



For this colour structure there are six singular regions (plus the ones generated from 
permutations like ki 



to be considered [|l^, [19 



first triple-coUinear region: kipi ~ k2Pi ~ A;i/c2 ^ 0; 
second triple-coUinear region: kip2 ~ k2P2 ^ kik2 —>■ 
double-soft region: ki,k2 ^ 0; 
first soft-coUinear region: fci — 0, k2Pi — >■ 0; 
second soft-collinear region: A;i ^ 0, k2P2 — ^ 0; 
double-coUinear region: kipi 0, k2P2 — ^ 0. 



In the first region the singularity is controlled by the coUinear factorization formula: 

\M,,^,,F{Pl,P2,k„k2,<P)\l ^ ^^''^'2'''^' ^^nt'L.fe)l-^SU(^3Pl,P2, 0)P , (76) 

where p^^^^ is the abelian part of the splitting function that controls the collinear 
decay of an initial state quark into a final state gluon pair. This function can be obtained 
from Eq. ( |A.6| ) with the replacement in Eq. (|59|). 

^^A factor 1/2 has been included to account for the two identical particles in the final state. 
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In the double-soft region the factorization formula is obtained by factorizing the two 
eikonal factors for independent gluon emissions (see Eq. (A. 3) of Ref. ||19||): 



\Mf^^^^APuP2,h,k2A)\L ^ (4vr/i2^«s)'l6C^5i2(A;i)5i2(A;2)|^;°U(0)r, (77) 
with Suik) defined in Eq. (H). 

In the soft-collinear region, say when k2 and kipi ^ we have instead (see 



Eq. (A.5) of PI): 



i-Mit..F(pi, P2, k,, k2, 0) r ^ (4v-as) %^ , , - 

Zi)piki{pi- ki)k2P2k2 
•|^SU((l-^i)pi,P2,0)P, (78) 

where Zi is the momentum fraction of the collinear gluon of momentum ki and can be 
identified with the one parametrizing the triple collinear splitting in Eq. ([TBI). Notice that, 
since the soft gluon of momentum k2 does not resolve the pair of collinear partons, there 
is no non-abelian contribution in Eq. (^^. 



In the double-coUinear region we have, when e.g. kipi —>■ and k2P2 0: 
\MfL,FiPuP2, ku k2, 0) P ^ i^T'n'^'-^ , (1-^1' ^) A. (1 - ^2, e) 



• I^SU ((1 - Zl)Pl, (1 - Z2)P2, 0) I' , (79) 

where zi and Z2 here represent the momentum fractions (see below) involved in the two 
collinear splittings. 



As it happens for the CpCa contribution, Eq. ( |76D supplemented with the definitions 
(|63|) is able to approximate correctly also the double-soft and soft-collinear regions in half 
of the phase space. But, at variance with the CfCa case, the same formula cannot describe 
correctly the double-coUinear region, since that one corresponds to the emission of gluons 
from different legs, i.e., with a kinematical configuration completely different from the 
triple-collinear case. Therefore, the strategy followed for the other colour factors does not 
work in this case. 

In order to overcome this problem there are in principle two strategies. The first one 
is to split the phase space in order to isolate the double-collinear region and perform the 
calculation separately for its contribution using the expression in Eq. (|79D. The second is 
to modify Eq. ([76|) in order to enforce the correct singular behaviour in all possible limits. 
We decided to follow the second strategy and for that we have first studied Eq. (^) with 
the definitions ( |63D and isolated the terms that do, incorrectly, contribute (terms with X2 
in the denominator) when the collinear gluons are emitted from the different legs. In this 
way, we were able to find a slight modification of Pgiglgs in Eq. ( |A.6| ) that allows to take 
into account the double-collinear region without spoiling the behaviour in the other regions 
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as 



f)(ab) ^ J ■"123 



X3 



1 + X3 xl+ X] 



+ 



S123 



X3(l -Xi) + (1 -Xs)^ 



' V,(^i)/,(^2)-e(l + e) 

-a;2 



Sis 



X1X2 



e(X]^ + X1X2 + X2 



X1X2 



f,iz2) + e\l + xs) 



Sl3 



+ (1^2) . 



^0) 



With respect to the expression of Pgiglqa of Eq. ( |A.6| ), the only difference is due to the 
introduction of the extra factors fq{z). The function fa{z), anticipating that a similar 
approach will be followed in the gluonic channel, is defined by 



faiz) 



-Paail-Z^e) 



(1^9 ■ 



^1) 



where Paa are the coUinear splitting kernels in Eqs. 



The function D^^^g^q^^ depends on the new momentum fractions zi and Z2 of the gluons 
with respect to the incoming antiquark of momentum p2- These momentum fractions should 
be the ones relevant for the double-coUinear limit. Our improved factorization formula is, 
outside the second triple-coUinear region given by 



\M 



(0) 
qq- 



■99 F 



{pi,P2,ki,k2,(j))\ 



ab 



where D^^^\ „ ,„ 
and by setting 



)(ab) 



ClD 



(ab) 

q^9i92{q-i) I 



|-MJ°U(^3Pi,P2,0)P 



^2) 



is obtained from Dg^g^q^, in Eq. (^) with the definitions in Eqs. (|5 



P1P2 



Z2 



P1P2 



^3) 



With Eq. (^) we can consistently approximate the relevant matrix element in the region 
from uq to Umax as we did in the other channels, keeping into account all the singular 
regions. In fact in the triple-collinear region zi,Z2 — * and fq{zi),fq{z2) 1. Therefore, 
in this limit Eq. ( ^2|) reduces to Eq. (|7^) . The factors fq{z) become relevant in the double- 



collinear region, since they ensure that the correct limit is recovered when piki 
P2^2 (and the same for ki ^ ^2). 



and 



Notice that the modification of the triple-collinear formula does not spoil the process 
independence of our calculation: it just allows to write an 'improved' formula that correctly 
interpolates all possible (double-) coUinear and soft singularities in the region of phase space 
where we have to integrate it. Therefore, with this approach we can avoid to split the phase 
space in regions where different approximations should be applied. 



It is worth noticing that the modification in Eq. ( PPD makes the calculation more involved 
already at the level of the angular integrals, mostly due to the introduction of the 'new' 
momentum fractions Zi and Z2- 
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For this colour structure we have to subtract the contribution from the factorization 
counterterm, which can be written as 



da 



FCT 



dq^dQ"^ 



as /" 2n ( da{pi,p2,(j),k) 

271 J ^^y''^^'F)\y dq^dQ-^dct) 



dq'^dQ'^ 



^4) 



where da{pi,p2, 0, k) corresponds to the cross section for the production of F and only one 
extra gluon (see Eq. pi)) and 



R(x,iJ,p) 



e^ii ^^^r(l-2e) 



^4 



2\ <^ 



i5) 



with P^g^{z 



n.(^'0) 



the regularized AP splitting function and fip 



the factorization scale. In the limit of small g-r and after taking moments with respect to 
z, the contribution from the counterterm factorizes as 



sSUt)(A') = 2Cj.^«('V,«) 



^6) 



where Tq^iN.e) and /C are defined in Eqs. (|35D and (|68D , respectively. Therefore, in the 
qT ^ ^ limit also the contribution of the factorization counterterm becomes process inde- 
pendent. 



Our final result for the 

(2) 



moment of the factorized contribution S 



W(ggab) 



{2)R 
ijg(ggab) 



(0) = 



(0)-S;;(^CT)(0) is: 



(0) =Cf}C 



^ + -)Cf J-,,-(0, e) + - / 2P,,(z, e) log z 



+ Ci 



- 21og' \ + 91og" ^ - 2 + -vr" log^ + 16C(3) - vr" - — 



^7) 



in agreement with the result of Ref. PD| for Drell-Yan. Notice that since 7gq''(0) = there is 
no contribution from the factorization counterterm to TiiJ, uJO). As we did for the other 



colour factors, we have isolated in the first line of Eq. 
by a similar term in the virtual contribution. 



_ 99(ggab) ^ 

'^) the part that will be cancelled 



A comment to the results obtained so far is in order. The formulae in Eqs. (|B7D, (|7DD, 
(fflp, ([75|), ( ^71) show that the contribution to Sqg(O) from double real emission are actually 
independent on the specific process in (|I|). This feature of the double real emission, which 
is due to the universality of soft and coUinear radiation, will persist also in the gluon 
channel. The explicit results obtained so far all agree with the ones obtained for Drell-Yan 



in Ref. El. 
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4.2 Virtual corrections 



The second part on the calculation of J]gg{N) involves the (one-loop) virtual corrections 
to single-gluon emission. The corresponding soft and collinear limits have been recently 



studied in Refs. pT| , |22| , |23| . The kinematics is the same as at 0{as) and the singularities 
originated by the same configurations discussed above Eq. (^). In the first collinear region 
the interference between the tree-level and one-loop contributions to single gluon emission 



behaves as 21, 22 



+ c.c. ~ 



ziPik 



X 



Pqq{Zl, e) [Mf^ip {ZlPl,P2, 4>) M'gg^^p {ZlPl,P2, 4>) + C.C, 



+ 2g 



2 ( 4vr/i2Y 



2pik J 



In Eq. (^8]) there are two terms. In the first one the tree- level splitting kernel Pqq{zi,e) is 
factorized with respect to the interference of the renormalized one-loop amplitude Ai^^^^p 



and the tree level one Ai 



(0) 

<?<?- 



The second term contains instead the unrenormalized one-loop correction to the splitting 



kernel Pg]l(g)„{zi,e) times the Born matrix element squared. The function P^^r ) {zi,e) 



(1) 



controls the one-loop collinear splitting of an initial state quark into a final state quark 
with momentum fraction zi, in the CDR scheme. Its explicit expression can be derived 



from the results of Ref. El], E2| and is up to 0{e 



,(1) / N Ct 



(An) 



1 1. . 1 



CAPqqix, e)[-- + - \og\l - x) + Li2 — Lis (1 - x) 



+CpPgg{x, e) ( -- log(x) - 2 log(x) log(l - x) + 2 Lia (1 - x) 



-CpiCp-CA)x 



^9) 



where 



r(i + 6)r^(i- 
r(i-26) 



(90) 



A factorization formula similar to Eq. 
state antiquark. 



holds when the gluon is radiated by the initial 



Let us now consider the soft region. At one-loop order, for a general amplitude with n 
hard partons, soft factorization formulae involve colour correlations between two and three 
hard momentum partons in the matrix element squared p3|. Nevertheless, in the case of 
only two hard partons the soft singularity is controlled by a simpler factorization formula 
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(see Eq. (57) of Ref. 
-^SLf (Pi,P2, K <P)M'~'^Z,FiPuP2, k, <P) + c.c. ^ 167ras/i'^ Cf 



Si2{k) (-MSip(0)-MSU(0) + C.C.) +S'S{k)\M^S-*M)\^ , 

(91) 



where 



1 r^(l-e)r3(l + e) 
e2 r2(l-2e)r(l + 2e) 



(92) 



is the unrenormalized one-loop correction to the tree-level eikonal factor. Likewise at 0{as) 
(see Eq. (^if)), in Eq. (|9ll) colour correlations are absent, and the factorization formula is 
similar in structure to the coUinear one in Eq. (p8|). 

One can verify that, as it happens at leading order, the factorization formula Eq. (^) 
with Zi = z = Q'^/s correctly reproduces also the behaviour in the soft-region, given by 
Eq. m. 



Furthermore, by expressing kpi in terms of qr and using Lorentz invariance as we did 
at leading order, a single factorization formula in the small qx limit is obtained: 



{pi,P2,k, 0) M'gjZgF {Pl,P2,k,(j)) +C.C. 



Anasfi^' 2{1 - z) 



X 



P,,i^, e) ( Mf^l^ (0) M^;^^ (0) + c.c. 



qq^F 
2" 



(93) 



and this formula can be used to approximate the virtual contribution in the full phase 
space. The same formula can be obtained by defining the coUinear momentum fraction in 
Eq. (PP as: 



kp2 
P1P2 



(94) 



and similarly when pi ^ P2- It is important to point out that, at variance with what 
happens in the double real emission contribution, here a process-dependent information 



appears, i.e. the one-loop matrix element Aigg^p {(f)). The most general structure of the 
product M^g^l^pAi^J^^p + c.c. is, according to Eq. (pSl): 



•'^'■qq^F-'^'-qq^F 



as [A'KLi'^Y rfl -e) 

+ C.C. = — 



27r V / r(l - 2e) 



ICf 3Cf 



+ \M 



,(0) I 



(95) 



In Eq. (^) the structure of the poles in e is universal and fixed by the flavour of the 
incoming partons, whereas, as discussed in Sect. ^ the finite part is parameterized by a 
scalar function ^^(0) depending on the kinematics of the final state particles. 
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The contribution from the UV counterterm in the MS scheme is: 



2 



(96) 



By approximating our matrix element with Eq. (|93|) , the calculation can now be performed 
quite easily as in Eq. (^). Using Eq. (pSj), and adding the contribution of the UV coun- 
terterm in Eq. ( p6D we find: 



E^f (0)=C^/C< 



\{2Cf + Ca)-- (ca log ^ + 3Cf + Po 



2 
e 



2p,g{z,e)\ogz + CA 



g2 39 



5 + -7r^ log^ + --27r^ + 2 log ^ - 3 X" (0) 



+ /3olog^ f3-21og^ 



(97) 



where is the renormalization scale at which as is now evaluated. The terms involving 
JF^q(0, e) and Pqq{z, e) in Eq. (|97|) are the ones that cancel against the corresponding terms 
in Eqs. (^), ( ffSD and (^). In the case of Drell-Yan, by using Eq. (p9D, our result agrees 
with the one of Ref. liC 



4.3 Total result for the qq channel 



After adding the real and virtual contributions in 



^^(0) — ^gg(nid) 

all divergent terms in e cancel out and we find: 

S(j(0) = log^^[-2C7|] 
Qt 

+ log" % [9C| + 2C^/5o] 

+ log^ 

(It 



(0) + yP-]^,. (0) + Yp-]!^ , (0) + s , , 



.(2)R 



(0) + S^'}^ . , 



(0) + e J^(0) (98) 



Cl ( -vr^ -71+ 2CMf (0) + CpCa ( ^ " y ) + ^^^/^^^ 



- 4C(3) ) - 'iCpA'M + C^fC^a ( - ^vr^ + 6C(3) 



'35 






3 


13 


11^ 


T ~ 


18 



9 



-TT 



(99) 



where we have set again [j?p = fi'j^ = Q^. It is worth noticing that the process dependence 
in Eq. (|9^) is fully contained in the function A^{(p). 
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Once one moment (the iV = in this case) has been computed, it is quite simple to 



extend the calculation to a general value of by studying the combination [jT3 



^{N) - S(0) = j dz (^^ - 1) 



dcr 



dao/d(f) dq^dQ'^ 



(100) 



Here, the factor (2;^ — l) eliminates singularities in the integrand when z — >■ 1 and allows 
to set qt = once the integral over the variable u has been done (in most of the cases it is 
possible to set qt = even before integrating over u). In that sense the complexity of the 
calculation is considerably reduced and the result can be expressed as 



E;|^(iV)=log3^[-2C|] 

Qrp 

■)2 



+ log" % [9C| + 2CfI3o - 6Cf7J;HA^)] 
It 



Qt 



7 



35 



TT 



CpCa I — — — 



+ 



18 3 

+ (2/3o + 12CF)7SHAr)-4(7«(iV)) 
15 \ / 13 11 



y 



1 



-4C(3) \^CpCa 



(Ar+l)(iV + 2) 
TT^ + 6C(3)^ - 3CMf (0) 

1 



^C^n^Tn (1 + ^vr^) + 27[^))(iV) + 2Cf^^^}[N) + 2 
+27(J)(iV)A^(<^)-2C^(/5o 



(iV + l)(iV + 2) 



3Cp 



(iV+l)(iV + 2) 2 



(101) 



where 7^^''^ ( A^) is the non-singlet space-like two-loop anomalous dimension |^ . The ex- 
traction of the resummation coefficients for the qq channel from Eq. ( [L01| ) will be performed, 
along with the corresponding one for the gg channel, in Section p. 



5 The calculation at 0(ag): the gluon channel 



The strategy for the computation of the 0(«s) contributions in the gluon channel is the 
same as the one developed for the qq case. In a similar way, we first consider the double 
real emission contribution and then the virtual correction. 

Let us first discuss the contribution coming from the factorization counterterm, that 
will be subtracted from the real corrections in the next subsection. By following the same 
steps that lead to Eq. (|^) we obtain 



^,WT)(^) = 2CA^..(iV,e) 
+ 2CF^,,(iV,e) 



e 

Ik 

e 



(102) 
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Eq. ( |102| ) contains two terms. The first one, due to the subtraction of one collinear gluon, 
is analogous to the one in Eq. ( p6D and contributes to both C\ and CaTr colour factors. 
The second term is due to the subtraction of a quark (antiquark) collinear to the initial 
state gluons, and contributes to the CfTr part. The function J^gq in Eq. ( |102| ) is defined 
as in Eq. (^Sj) by 

C,Tg,{N^e) ^ n^"" dzz^^^Z$M4^ . (103) 



^{l-zY-Azqi/Q 
The — >• limit can be safely taken and Tgq{N, e) gives 

^,,(iV, 6) - 2 dz z^ (^ii±ii^ _ ez^ = 2^a)(iV) - 26^ . (104) 



5.1 Real corrections 

The contributions to H^^g from double real emission fall in two classes: 



g + g g + g + F . 



The kinematics is the same as discussed at the beginning of Sec. |4.1| . As we did for the 
quark channel, we will first perform the calculation for a fixed moment and then extend 
it for general N. Since the N = moment is divergent for the gluon channel (see e.g. 
Eq ( |104| )), we start from = 1. Furthermore, as it happens at LO, spin-correlations 
appear in the collinear decay of a gluon. Nevertheless, since the correlations cancel out after 
integration, we will use in the collinear factorization formulae directly the spin-averaged 
splitting functions. 



5.1.1 Contribution from qq emission 

For this contribution the strategy followed for the CpT^ and CfCa terms in the quark 
channel applies. The singular regions are: 



first triple-coUinear region: kipi ~ k2Pi ^ kik2 0; 
second triple-collinear region: kip2 ~ k2P2 ~ ^1^2 —>■ 0; 
double-soft region: ki,k2 



In the first triple-collinear region the factorization formula reads 

\Mf^^qgAPuP2,k„k2A)\' ^ ^^''^'^"'^' ^.^.-...(.3)|A^SlF(^3Pi,P2,0)r, (105) 
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where Pg^g^^g^^g^) is the sphtting function that controls the decay of an initial state gluon 
into a final state quark-antiquark pair and a gluon. It can be obtained from the expression 
of Pqiq2gz that dcscrlbcs the decay of an (off shell) gluon into a final state quark antiquark 
pair plus a gluon, given in Eq. ( |A.8| ), with the crossing transformation ([59|). In the double- 
soft region the factorization formula is the same as in Eq. ( |6TD with Cp —* Ca and, likewise 
in the quark channel, the soft behaviour is correctly taken into account by Eq. ( |105|) with 
the definitions (|63|). Therefore, we can follow the strategy successfully applied in the quark 
channel to obtain 



t^r- = CpUfTn ;C f- A _ 1^ log ^ + 1^ log 91) 



(106) 



and 



SL(qqnab)(l) = CAUfTn /C (^---^,,(1, e) + - - — + - log ^ - - log ^ 

4, 11 .1 



In Eqs. ( |106|jl07| ) we have already subtracted the CpTji and CaTr terms from the factor- 
ization counterterm (with = 1) in Eq. ( |102| ). Furthermore, in Eq. ( |107| ) we have isolated 



a divergent term that will be cancelled by a similar term in the virtual contribution. The 
explicit expression of the function Tgg{l,e), defined in Eq. (0), reads 



^,,(l,e) = -^ + 21og^. (108) 

6 Qrp 



5.1.2 Contribution from gg emission 



The calculation of the gg C\ contribution to Tjfg parallels the one for the part in 
the quark channel since the singular configurations have the same complicated pattern 
as described above Eq. ([T6|). The triple-coUinear region is controlled by the factorization 
formula 

\Mfg^gg,{p,,P2MMA)? ^ ^^^^^^^^^^^^ [m 

where the function Pg^g^g^(^g.j^) is now obtained by applying the crossing transformation (^) 
to the splitting function Pgig2g3 that controls the coUinear decay of a gluon into three final 
state gluons, given in Eq. ( |A.11| ). 

The factorization formulae in the soft-collinear and double-collinear regions are analo- 
gous to Eqs. ( [78| , [r9[ ), and can be obtained from them by conveniently changing the colour 
factors {Cp Ca) and splitting functions {Pgg Pgg)- The factorization formula in the 
double-soft region receives two contributions analogous to the ones in Eqs. ( [73|j7^ ). 

Likewise in the quark channel, Eq. (|109|) with the momentum fractions defined as in 
Eq. (^) approximates correctly, again in half of the phase space, all possible infrared 
configurations but the double-collinear one. 
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In order to proceed further, we use the technique developed for the quark channel. As 
before, we first study the behaviour of Eq. ( |109| ) with the definitions (p3| ) in the double 
coUinear limit and identify the terms that do (incorrectly) contribute in that limit. Those 
terms have to be modified in order to enforce the correct double-collinear limit without 
affecting the singular behaviour in the other regions. The modified splitting function we 
obtain is: 



n _ pnon-sing , 

-^919293 919293 ~'~ 



psing-1 I / psing-1 , psing-l^ /■ 

919233 ' \ 919392 ~ 1:^(11 ci-2 I Jgy'^i 



939192 



+ (1^2) 



(110) 



The first term 



^non— sing ^2 



919293 



+ 



'123 



:i-e) 



4,3 + ^(1-^) + ^ 



■Sl2'Sl3 



, X1X2 — 1 3 5 

4-!-^ 1 h -X3 

1 - X3 2 2^ 



^ , xi(l + 2xi) l + 2a;i(l + xi) 
X2X3 - 2 H h 



2(1-X2)(1-X3) 



+ (5 permutations) 



111) 



contains the part of Pgig2g3 in Eq- ( |A.11|) that does not contribute to the double-collinear 
limit. Therefore, this part of the splitting function does not need any modifications. The 
variable tij^^ is defined in Eq. ( |A.4| ). 

The second part is instead modified with the introduction of the function fg{z), defined 
in Eq. (|HT|). The functions P^j^f^g^ and Pg^g^'^^ are 



psing-l _ ^2 £123 / 3^13^2 - 2 (1 - X3(l - ^3))^ 
-^919293 '-^4 1 



S12 



X3 



X3Xi(l - Xi) 



(112) 



psmg-2 _ ^2 ^123 f ( 1 - X2) ( 1 - 2X3) l-2Xi(l-Xi) » 

™^ ^S,2SrA X3(l-X3) + 2X2X3 ^ 



Our improved factorization formula is therefore 

\Mf^^^^APuP2,kuhA)\' ^ ^^^^^^^^D^^^^ (114) 

As for the quark channel, the expression of IDg-,g^g2{^g^) is obtained from the one of Dg^g^g^ 
in Eq. (|11CI|) by using Eqs. (^91), (|63D and defining zi and Z2 through Eq. (|83[) . 

In the triple-coUinear limit fg{zi), fg{z2) —>■ 1 and the various contributions in Eq. ( |110| ) 
reconstruct the triple-collinear splitting function -Pgig293- The role of the functions fg is 
again to enforce the correct behaviour in the double-collinear region. It is worth stressing 
that there are in principle many ways to conveniently modify the splitting function and 
that we have tried to find the simplest one that fulfills all the requirements and can be 
integrated afterwards. 
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The function Pg^g^g^ in Eq. ( |A.11|) has by itself the most comphcated expression among 



the various Paia2az because one has to sum over six permutations. Besides that, the mod- 



ification in ( p.lO|) makes the angular integration very involved. Since many of the ensuing 
terms have an additional singularity as S2 ^ some of the angular integrals in Ref. |^ 



have to be evaluated one order higher in e. Once the angular integrals have been performed, 
one has to face an additional complication: 'spurious' l/^'l^ and singularities appear 
in the intermediate steps, which of course cancel in the final result, but create additional 
problems to take the g-p limit. The final (factorized) result is 



gaiss) 



;i) = Ca/c 



3 311 1, 
- + -— + - log — 



e 6 



99 
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2zPgg{z,e) hgz 



log — + ^ log — 



82 



An' 



llvr^ 



10C(3) 
(115) 



where we have isolated in the first line the terms that will be cancelled by analogous virtual 
contributions. Notice that, since 7gg (1) = — |n/TR, there is no contribution to Eq. ( |115| ) 
from the factorization counterterm in Eq. ( |102|) . 



5.2 Virtual corrections 



(2) 

We finally compute the small-g^ behaviour of the virtual contribution to S^g . The calcu- 



lation parallels the one for the quark in Sec. ^4.2| , and the singular configurations are the 
same as at leading order. For the coUinear limit, say when kpi — > 0, we can write a formula 
similar to Eq. (Rl) 



M^^^l„p (pi,P2, fc, 0) M)!^JZnF (Pl.P2, k, 0) + C.C. 



{l)u 
99-^9 i 



2e 



X 



((0)t 
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Pggizue) iM\'lp {ZlPl,P2, 



.2 /'47r/i V 5(1) 



2pik J 



zipik 

gg^F 

f(^iP1'P2, 



M^}}^„ {ZiPi,P2, 0) + C.C. 



99 



(116) 



and pj^^^^^glzi, e) is the un- 



where Pgg{zi,e) is the tree-level splitting kernel in Eq. ([ 

renormalized one-loop correction to the AP kernel for the coUinear splitting of an initial 
state gluon into a final state gluon with momentum fraction zi, in the CDR scheme. Its 
explicit expression can be derived from the results of Ref. |^l| and is up to 0{e^): 



9^(9)9^ ' ' 



:i-x 



{Any 



12 vr' 
^99(3;, e) ( — 2 log(3;) - 2 log(l - x) log(x) + — 



--{CA-2nfTR)x 



;il7) 



A similar formula holds when the gluon is radiated by the initial state antiquark. In 
the soft region, the factorization formula is the same as in Eq. (|9TD with Cp Ca and one 
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can verify that, as it happens in the quark channel, Eq. ( 116 ) with Zi = z = Q'^/s correctly 
reproduces also the behaviour in the soft-region pT| , |23| . In the same way as for the quark 
channel we can write down a single factorization formula in the small qx limit: 

Mf^l,, {PUP2, k, k, </>) + c.c. ^ 4vras/x^-2(l-.) 

Qt ^ 

(Ks^F (0) -M^^l^ (0) + C.C. 



X 
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/ 47r/i' 
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2\ ^ 



zYP 
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9^(9)9 



118) 



According to Eq. ( pH]) the renormalized amplitude A4^gg_^p can be written, up to C(e°) as: 



w(0)t , _ ^ 



as /47r/i2Y r(l-e) 
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^+<(0))|A^ 



(0) I 
99^F\ 



(119) 



In the case of Higgs production, in the mn -C rritop limit and including also the finite 
renormalization to the effective ggH vertex , the function A^{4>) is given in Eq. 



The contribution from the UV counterterm (in the MS scheme) needed to renormalize 



the splitting kernel P 



(1) 

9-*{9)9 



is0: 



-(2) 
^ggUVCT 



{N)=CATgg{N,e)lC 



--)/5o. 



(120) 



By approximating our matrix element with Eq. ( 118 ), using Eq. (|119D , and adding the 
contribution from Eq. (|120|) we find 



S(f (1) = Ca/C< 



+ CA{-\oi% 



e 

2 + 



11 



Qt ^ 
65 Svr^ 
18 ^ ~ 



-^^99(1'^) - - / 2 zPgg{z,e) log z 



Q2 11^2 



+ ^.,T, + <(,)(-H + 21og|)-/.olog|(-^ + 21og|)}. (121) 

The terms involving Tgg{l, e) and Pgg{z, e) in Eq. (|121|) cancel the corresponding divergent 
contributions in Eqs. ( |107|) and ( |115| ). 



5.3 Total result for the gluon channel 

After adding all the contributions in 

^^W=^SW^+ ^S:.na.)(l) + S(^^)(l) + S(f (1) , (122) 

^■^The total UV counterterm in the case of Higgs pro duction would be three times this one, but we have 
included part of it in the renormalized amplitude ( |119|) . 
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all divergent terms cancel out and we obtain 



E(J(l) = log»^[-2Ca 
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+ log2^[8CA/3o + 4CAn/TR] 
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+ log — 
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229 TT^ 



16 



- ) + -CAUfTn - -CpUfTn + 2CaA^ (0) 
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Cl 2C(3) 



16 11 



-TT 
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(123) 



where we have set again = /x^ = 



The contribution for general can be computed as explained in the previous section 
for the quark channel. The total result is: 



HT 



+ log2^[8CA/5o-6CA7SHiV)] 
Qt 

+ log^ 
(It 
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C\ (f + y ) - ^^CAU^Tn + 2CaA^{<P) 
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+27f)(iV) + 2V;HiV) ( ^f(0) + Ca- ) +4C^n;7^^^(iV) 
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(124) 



Here 7gg (iV) is the singlet space-like (gluon-gluon) two-loop anomalous dimension [12 
whereas the coefficient rjAi^ has origin on the N moments of —Pg^i^z), see Eq. ( |104| ). 



6 Final results and discussion 



We can now compare the results obtained in the previous sections with the second order 
expansion of the resummation formula in Eq. ([T8|). As for the A^-dependent contributions 
in Eqs. ( |101D , ( |124| ), they fully agree with the ones in Eq. (|15)Q. This agreement can be 
considered as a non-trivial check of the validity of the resummation formalism, because the 
expressions in Eqs. ( |101[ ),( P^2^ are completely general and the process dependence is fully 
embodied in the functions A^{(j)). As an alternative, given the resummation formalism for 
granted, the result in Eqs. ( |TOl]) , (|124|) can be considered as an independent re-evaluation 
of the two-loop anomalous dimensions. 
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We have checked that the results in the quark singlet channel are also in agreement with Eq. (Jig). 
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As far as the A^-independent part is concerned, it can be used to fix the coefficients A^"^^ 
and B^^\ By comparing the single-logarithmic contributions in Eqs. ( |101[ ), ( |124| ) with the 
one in Eq. (|T8|) we obtain for the coefficient A^a^: 

Ai2)=irA« a = q,g, (125) 

where K is given in Eq. (^TJ), thus confirming the results first obtained in Ref. |Ty, |13[. By 
comparing the non logarithmic terms we find that the coefficient B^'^^ can be expressed as 
well by a single formula for both channels: 

^ = -2 7f +/5o (^C^aTT^ + ^r(0)) a = q,g, (126) 

where ^y!^^ are the coefficients of the 6{l — z) term in the two-loop splitting functions Paa\z) 
4l| , given by 



c|(^^.aa3,).c.c.(lI.Hf!-3a3,)-c.„,r„(i.|!) 

Cl + 3C(3)) - CpnjTn - ^CAUfTn . (127) 



From Eq. (|126D we see that 5*^^^ besides the term which matches the expectation 

from the 0{as) result, receives a ^rocess-c?ej>en(ieni contribution controlled by the one-loop 
correction to the LO amplitude (see Eq. (|38|) ). Thus, as anticipated at the beginning, 
although the Sudakov form factor in Eq. (|^) is usually considered universal we find that it 
is actually pro cess- dependent beyond next-to-leading logarithmic accuracy. 

However, by using the general expression in Eq. ( |126| ) it is possible to obtain B^"^^ for 
a given process just by computing the one- loop correction to the LO amplitude for that 
process. For the Drell-Yan case, by using Eq. (^), our result for T.qq (N) agrees with the 
one of Ref. [0, confirming the coefficient Bf^^^ in Eq. (0). 

In the interesting case of Higgs production in the rxitop oo limit, by using Eq. (^) 
we find[3 : 

^ = ^1 (f + f - 6C(3)) + ^C, nj Tn - uj T« Q + ^tt^) - . 

(128) 

In particular, this result allows to improve the present accuracy of the matching between 
resummed predictions |43| and fixed order calculations [J5]. 



(2)H 

Notice that in this case, the coefficient Bg turns out to be numerically large. Actually, 
for Uf = 5 we have Bg /Bg ^ —14, whereas for Drell-Yan the same ratio leads to 
^^2)Dy^^(i) ^ —1.9^ i.e. about 7 times smaller than for Higgs production. Both the 
appearance of a C\ term (compared to Cp in the quark case) and the size of the one- 
loop corrections to Higgs production are the reasons for the large coefficient. Clearly, the 

^^Actually, using the results of Ref. Q for the two-loop gg ^ H amplitude, one can also obtain Bg'^^^ 
for arbitrary rritop- 
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use of Bg^^^ in the implementation of the resummation formula will have an important 



phenomenological impact Actually, one can expect that the inclusion of Bg^^^ , which 
will tend to reduce the resummed cross section, will partially compensate the increase in 
the normalization produced by the (also) large coefficient Cg^J^ 

The fact that the Sudakov form factor is process-dependent is certainly unpleasant. 
Usually it is called the quark or gluon form factor, since it should be determined by the 
universal properties of soft and collinear emission. With the result in Eq. ( |126| ), instead 
we find, for example, that the form factor for 77 production is different from the one for 
Drell-Yan. Moreover, since the hard function depends in general on the details of the 
kinematics of F (in case of 77 production it would depend, e.g., on the rapidities of the 
photons), the same happens to the coefficient B^'^^ and thus to the Sudakov form factor in 
Eq. (D. 



However the results in Eq. (|^) and Eq. ( |126| ) suggest a simple interpretation We 
can see in Eq. (^) that the pro cess- dependent coefficients functions C^'^^{z) have two 
contributions. The first has a collinear origin and is driven by the 0{e) part of the Pab{z, e) 
kernel (see Eqs. (^7|)). The second has instead a hard origin, and contains the finite part 
of the one-loop correction to the leading order subprocess. As a consequence, the scale at 
which as should be evaluated is different for these two terms. In the collinear contribution 
as should be evaluated at same scale as the parton distributions are, i.e. bl/b"^. By contrast, 
the correct scale at which as should be evaluated in the hard contribution is the hard scale 



As discussed in Ref. this mismatch, that affects the resummation formula in its 
usual form Eq. (H), can be solved by introducing a new pro cess- dependent hard function 
H^[as{Q'^))- The ensuing resummation formula is 



W[,{s- Q, 6, 0) = V / dz, [ dz2 Cca{as{bl/b''),z,) Ca(as(&^/6'), Z2) S{Q^ - z^z^s) 
„ Jo 



daUQ',<ys{Q' 



Sc{Q,b) 



(129) 



where 



da^Q^asiQ' 



da, 



{LO)F 



(130) 



As discussed in Ref. p^, this modification is sufficient to make the Sudakov form factor 
Sc{Q,b) and the coefficient functions Cabias^bl/b"^), z) process-independent, with Cab and 
being dependent on the introduced 'resummation-scheme'. We point out that this 
modification is not only a formal improvement, since, once a resummation scheme is fixed. 



the resummation coefficients in Eq. ( |129[ ) are now universal and it is enough to compute 
the function at the desired order for the process under consideration. 

Summarizing, in this paper we have exploited the current knowledge on the infrared be- 
haviour of tree-level and one- loop QCD amplitudes at 0(ag) to compute the logarithmically- 
enhanced contributions up to next-to-next-to-leading logarithmic accuracy, in an general 
way, for both quark and gluon channels. Comparing our results with the gr-resummation 
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formula we have extracted the coefficients that control the resummation of the large loga- 
rithmic contributions. We have presented a result that allows to compute the resummation 
coefficient B^"^^^ for any process, by simply knowing the one-loop (virtual) corrections to 
the lowest order result. In particular we have obtained the result for the case of Higgs 
production in the large rritop approximation, which turns out to be numerically relevant for 
phenomenological analyses. 

The results of our calculation clearly show that the Sudakov form factor is actually 
process dependent within the conventional resummation approach. An improved version 



of the resummation formula where this problem is absent has been presented in Ref. [27 
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Appendix: Triple-collinear splitting functions 

In this Appendix we collect the various expressions of the triple-collinear splitting functions. 
Denoting by ri, r2 and the momenta of the final state partons that become coUinear, 
the triple-collinear splitting functions depend on the invariants Sij = (r^ + r^)^, S123 = 
■S12+S13+S23 that parameterize how the coUinear limit is approached, and on the momentum 
fractions Xj {i = 1,2,3) involved in the coUinear splitting. The splitting function for the 
coUinear decay of a quark q in qq pair plus a quark is 
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and the variable tij^k is defined as 
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The sphtting function for the q — > qgg decay can be decomposed according to the different 
colour coefficients: 
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and the abelian and non-abelian contributions are 
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When a gluon decays collinearly, spin-correlations are present. Here we are concerned only 
with spin-averaged splitting functions. When the gluon decays in a qq pair plus a gluon 
the splitting function is 
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In the case of a gluon decaying into three colhnear gluons we have: 
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